Modifying quantum walks: A scattering theory 

approach 



We show how to construct discrete-time quantum walks on directed, 
Eulerian graphs. These graphs have tails on which the particle making 
the walk propagates freely, and this makes it possible to analyze the walks 
in terms of scattering theory. The probability of entering a graph from 
one tail and leaving from another can be found from the scattering matrix 
of the graph. We show how the scattering matrix of a graph that is an 
automorphic image of the original is related to the scattering matrix of 
the original graph, and we show how the scattering matrix of the reverse 
graph is related to that of the original graph. Modifications of graphs and 
the effects of these modifications are then considered. In particular we 
show how the scattering matrix of a graph is changed if we remove two 
tails and replace them with an edge or cut an edge and add two tails. 
This allows us to combine graphs, that is if we connect two graphs we 
can construct the scattering matrix of the combined graph from those 
of its parts. Finally, using these techniques, we show how two graphs 
can be compared by constructing a larger larger graph in which the two 
original graphs are in parallel, and performing a quantum walk on the 
larger graph. This is a kind of quantum walk interferometry, 

1 Introduction 

Quantum walks are quantum versions of random walks. In both, a particle on a 
graph moves through the graph as time progresses. In a classical random walk, 
the path that the particle takes at a given time is determined by probabilities, 
while in a quantum walk it is governed by probability amplitudes. The result 
is that in a classical random walk, the motion is diffusive, while in a quantum 
walk, the motion is more akin to wave propagation. These walks were first 
proposed and studied by Aharanov, Davidovich, and Zagury jT]. They were later 
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rediscovered by a number or workers who were interested in them as possible 
sources of quantum algorithms [2 [3l |4] . The search for walk-based algorithms 
has been successful, and search algorithms 0, subset-finding algorithms [HI [7], 
a quantum algorithm that can solve a particular oracle problem exponentially 
faster than is possible with any classical algorithm [5], and, most recently, a 
quantum algorithm for evaluating the NAND tree [51 Unji have been found. 
There has now been considerable work on the properties of quantum walks, and 
some of it is summarized in two relatively recent reviews [111 I12j . 

Quantum walks come in two varieties, discrete and continuous. Here we 
shall consider only discrete walks. In these walks, there is a unitary operator 
that advances the walk one time step. In most versions of discrete-time walks, 
the particle making the walk is located on the the vertices of the graph, and 
states corresponding to the particle being located at a particular vertex form an 
orthonormal basis for a Hilbert space Hv, whose states describe the location of 
the particle. In order to guarantee the unitarity of the time-step operator, it is 
necessary to enlarge the Hilbert space by adding a quantum coin. For example, 
if the walk is taking place on the line, the coin space, Tic is two dimensional. It 
is spanned by the orthonormal basis and if the coin is in the state 

\R) the particle moves to the right on its next step, and if the coin is in the state 
\L) it moves to the left. On a more complicated regular graph, the dimension 
of the coin space is larger, and it is spanned by an orthonormal basis, each of 
whose elements corresponds to a direction. The quantum walk takes place on 
the space TL^^TLc- 

Here we shall consider a discrete-time quantum walk in which the particle 
is located on the directed edges, rather than the vertices, of the graph. The 
properties of a class of walks of this type, in which for each directed edge going 
from a vertex vi to vertex V2, there is a corresponding edge going from V2 to 

were explored in [13j and [14j . They have the advantage that a coin space 
is unnecessary, and that it is simple to define them for any graph. In this case, 
the Hilbert space that describes the walk is spanned by an orthonormal basis 
whose elements correspond to directed edges. That is, there are two orthogonal 
states corresponding to to each edge; one corresponding to the particle being 
on the edge going in one direction, and the other to the particle being on the 
same edge but going in the opposite direction. In |14j we considered walks of 
this type on a general graph, G, connected to two tails. Each tail is a half 
line consisting of an infinite number of edges, with one end going off to infinity 
and the other attached to a vertex of G. The particle propagates freely on the 
tails, for example, if it is on one edge moving to the right at one step, after the 
next step it is on the edge to the right of the one it was on and still moving 
to the right. The motion of the particle in the graph G is more complicated. 
This arrangement allowed us to study quantum walks from the point of view 
of scattering theory. Scattering theory was first applied to quantum walks by 
Farhi and Gutman, in the case of continuous-time quantum walks [5]; it our 
case it is applied to discrete-time walks. A freely moving particle approaches 
G on one tail, scatters in G, and has some amplitude to be reflected from G 
back onto the tail from which it came and another amplitude to be transmitted 
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through G onto the other taiL The properties of a walk in which a particle 
starts on one tail and is later measured to be on the other one can be found 
from a transmission function that is characteristic of the graph G. There is a 
corresponding reflection function that describes walks that begin and end on 
the same tail. Both are functions of a complex variable and are analytic in a 
region including the unit disc. 

Here wc would like to extend that work. We shall flrst consider graphs with 
directed edges in which there is not necessarily a directed edge from to vi if 
there is one from v\ to V2- The graphs will, however be Eulerian, that is each 
vertex will have the same number of edges entering it as leaving it. In addition, 
we shall allow an arbitrary number of tails. This leads to a transmission matrix 
instead of a transmission function. We shall then consider how the transmission 
matrix of the graph changes when the underlying graph changes. In particular, 
we shall see what happens when two tails are discarded and replaced by an 
edge connecting the two vertices to which they were attached or an edge is cut 
and replaced by two tails. The transmission matrix of the new graph can be 
calculated from the transmission matrix of the original graph. We shall also be 
able to combine graphs. In particular, if we have two graphs with tails, we can 
remove two tails, one from each graph, and replace them by a single edge that 
connects the vertices to which they were attached, thereby connecting the two 
graphs. The transmission matrix for the combined graph can be expressed in 
terms of the transmission matrices of the two original graphs. Finally, we shall 
show how two graphs can be compared by constructing a larger graph from 
them in which the two original graphs are in parallel. This allows us to do a 
kind of intcrferomctry on graphs. 

This approach has the advantage that it allows us to construct quantum 
walks on larger graphs from walks on smaller ones. A quantum walk is charac- 
terized by the transmission matrix of that graph. What we show how to do is 
compute the transmission matrix of a larger graph from those of smaller graphs 
that are its constituents. 

2 Basic formalism 

We shall begin by defining a graph in a rather general way. A graph G consists 

of a set V of vertices and a set E of directed edges, and two maps, i : E ^ V 
and t : E ^ V. These maps associate to each edge, e a point i{e), which we 
shall call the initial point, and a point t{e), which we shall call the terminal 
or end point. Wc allow both loop edges in which i{e) = t{e), and multiple 
edges, that is distinct edges with the same initial and terminal points. This 
type of ensemble is often called a digraph in the literature; all our graphs will 
be digraphs so wc will drop the "di". This abstract dcflnition has the obvious 
geometric realization in which we first embed the vertices as points in Euclidean 
three space, and embed the edges, each of which is a distinct copy e = [Og, le] 
of the imit interval directed from Og to Ig by mapping Og onto i{e) and Ig onto 
t{e). In the geometric realization an edge becomes a path joining i(e) to t{e) and 
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inherits its orientation from the unit interval. Thus loops carry an unambiguous 
orientation. For each v G V let uiy = {e\t{e) = v} and Ty = {e|i(e) = v}. These 
are, respectively, the sets of incoming and outgoing edges at v. Note that the 
sets and are disjoint for vi V2, as are the sets t^^ and Ty^. We also 
have that E = [j^^y ujy = [j^^v '^v 

We shall be interested in graphs that satisfy the condition = |t„|. 
Graphs with this property are called Eulcrian. If a graph is to be the underlying 
graph for the quantum walks we wish to study, it must be Eulerian. 

The basic picture of our quantum walk is the following. The particle making 
the walk is located on the edges of the graph, and when it passes through a 
vertex it scatters. If the particle is on an edge between the vertices vi and V2, 
it can either be going from vi to V2 (corresponding to the directed edge with 
i(e) = Vi and t{e) = U2), or it can be going from V2 to vi (corresponding to 
the directed edge with i{e) = V2 and t{e) = vi). If it is going from vi to V2, 
the next time step will carry it through V2 onto one of the edges leaving V2 ■ In 
order to define the walk we need a Hilbert space that describes a particle on 
the directed edges of the graph, and a unitary operator that advances the walk 
one time step. This operator is constructed from operators that describe the 
scattering at the individual vertices. 

We first construct the Hilbert space for the quantum states of a particle 
moving on the graph. Let Cly and Ty be the Hilbert spaces generated by taking 
the elements of LUy and Ty, respectively, as orthonormal basis elements. Let 
Uy : fly ^ Ty bc thc local scattering operator, and we assume that Uy is 
an isometry. By combining these local operators we are able to construct a 
unitary operator that advances the quantum walk one step. In particular, we 
define U : L^{E) L^{E), where L^{E) = ^^^yUy = ^y^yTy, such that 
U\fly = Uy. We call such a unitary U a quantum structure on the Eulerian 
graph G. 

Let G = {V, E, i, t) be a graph. The reverse graph G_r = {Vr, Eh, ir, ta) is a 
graph where the set Vr = ^ and Er = E, while ifl(e) = t{e) and tfl(e) = i{e). 
In this new graph Gn we have reversed the orientation on all the edges of the 
the geometric realization of G. We see that (r_R)„ = uiy and (wii)„ = r„. If the 
original graph G is Eulerian then so is Gr . If J7 is a quantum structure on 
G we will now define a reverse quantum structure Ur on Gr. Let L^(Er) be 
the Hilbert space generated by taking the oriented edges of Gr as orthonormal 
basis elements. Define the conjugate linear map R : L^{Er) — > L^{E) so that 
if \v2, vi) € L^{Er) is the basis element corresponding to the edge e considered 
as an edge of Gr, so that ^^(e) = V2 and tR{e) = Vi, then 



where \vi,V2) is the basis element corresponding to e considered as an edge of 
G, so that i{e) — vi and t{e) = V2- R is the orientation reversing map, and it 
is an isometry from {flR)v onto Ty and from (T/j)^ onto Qy, for all v gV. The 
reverse quantum structure on Gr is defined by 



R\v2,Vl) = \Vi,V2), 



(1) 



{UR)y = R-^U-^R, 



(2) 
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which is a linear isometry from {^Ir)v onto {Tj^)v The operator R defined here 
is closely related to the time-reversal operator for quantum walks defined in Ref. 

A graph G = {V,E,i,t) possesses a pairing if there exists a fixed-point- 
free involution A on the edges E such that t{Ae) = i(e) and i{Ae) = t{e). The 
initial and terminal points of e are the terminal and initial points of Ae. A graph 
possessing a pairing is clearly Eulerian. The simplest case of such a graph is 
when each pair of vertices which are connected at all have exactly two edges 
joining them one in each direction, a divided highway. Such graphs are said to 
be simple. Quantum walks on simple graphs were treated in [14) . 

LetG = iV,E,i,t) andC = {V\E\i\t') be a pair of graphs. Let = (/, F) 
be a pair of maps, f : V V and F : E ^ E' . </> is a graph morphism if 
i' o F — f o i and t' o F — f o t. If / and F are bijections we call a graph 
isomorphism. If is a graph morphism then F : ujy ^ ^f(v) and F : '''/(f)- 
F thereby extends to a linear map also denoted by F which maps into ^f(v) 
and also maps T„ into Tf(v)- Clearly F : L'^{E) L^{E') and < 1. Let G 
have a quantum structure U and let G' have a quantum structure U' . A graph 
morphism (j) is said to be a quantum graph morphism if 

FoU ^U' oF. (3) 

If (/) is a graph isomorphism or a graph automorphism which commutes with 
quantum structures as above then we call them quantum isomorphisms or quan- 
tum automorphisms respectively. 

In order to apply scattering theory to our walk, we need regions where the 
particle propagates freely, and no scattering takes place. For this reason, we 
will attach semi-infinite lines, or tails, to our graph. Let G be a graph with a 
finite number of edges and vertices, and Tic be the Hilbert space spanned by 
the states corresponding to its directed edges. We shall single out two subsets 
of the vertices, {vk\k — 1,...K}, and {ui\l — 1,...L} where we will attach 
incoming tails and outgoing tails Yi respectively. We shall make no further 
assumptions on these vertices. A given vertex may appear many times in one 
or both lists. We want to be able to clearly identify the distinct incoming 
and outgoing tails. The vertices of of the incoming tail Xj are denoted by 
kj, where k = 1,2, . . ., and the directed edges are \kj, {k — = \k,k — l)j 
for k > 2 and the attaching edge is \lj,Vj) — |l,0)j. The vertices of the 
outgoing tail Yi are denoted by mi, where m = 1,2,... . The oriented edges 
are |m;, {m + l)i) ~ |m, m + 1);, and the attaching edge is \ui, li) = |0, l)i. In 
this way given G and the two sets of vertices we can construct a new graph 
r = (G, (wi, ...,vk){ui, ...,ul)) where the vertices of new graph are those of G, 
the vertices kj 1 < k < oo for each tail X^, and the vertices {m)i 1 < m < oo 
for each outgoing tail Yi. The edges of F are the edges of G and the edges 
\k, k — and |m, m + l)i as above. In order that F be the underlying graph 
for a quantum structure it must be Eulerian , which implies that K — L, i.e. 
the number of incoming tails is the same as the number of outgoing tails. We 
call such a graph F an Eulerian graph with tails. 
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If wc arc given an Eulcrian graph with tails T, we want to study quantum 
structures [/ on F with the additional property that 

U\k + l,k)j = \k,k-l)j (4) 

on each edge of an incoming tail, and 

U\m-l,m)i = \m,m+l)i (5) 

on each edge of an outgoing tail. We say such a quantum structure is free. 
The "particle" freely propagates towards G along a incoming edge and freely 
propagates away from G along an outgoing edge. 

We now want to consider eigcnstates of U that correspond to the following 
situation. A particle approaches G on the tail X)., scatters in G, and then has 
amplitudes to leave G on any of outgoing tails Yi . A quantum state of this 
form is given by 

OO K QO 

;=o 1=1 m=0 

and satisfies the equation 

U\MO)) = e-''\MO))- (7) 

The first part of \tpk(d)) corresponds to the incoming particle, the second to part 
of the state inside G, and the final part to the outgoing particle. As we shall 
later show, the functions t^i'\0) and |V'G,fe(^)) are restrictions to the unit circle 
of functions that are analytic for |2| < 1 + e, for some e > 0, while IV'fe(^)) itself 
is the restriction to the unit circle of an analytic function from the punctured 
disc, < 1^1 < 1 + e into L°^{E). 

In order to begin the extension of the above eigenstate into the complex 
plane, define 



=0 



.{z)) = Vz-'|/ + l,Oi. (8) 



Note that 



U\aj+{z)) = l(|cT,.+ (^))-|0,l),), 

U\aj_{z)} = i|a,_(z)) + C/|l,0),-. (9) 



We then define 



K 

\Mz)) = Wk-iz)) + li^aA^)) +J2tf\z)\aj+{z)), (10) 
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such that it is the solution to the equation 



zU\M^)) = l^fc(z)). (11) 
This equation wiU be satisfied if and only if 

[/(|^G,fcW) + |l,0)fc) = -J.HG,k{^^)) 

+ (z)|0,l),). (12) 

This will be our key equation, and we shall now analyze it in more detail. 
Let Pg be the orthogonal projection onto Hg, and set 

\wk) ^ PGU\l,0)k- (13) 
If we now apply Pg to both sides of Eq. (fT2l) , we have that 

PgUI^gA^)) + \wk) - ^I^gAz))- (14) 
Defining Ug — PgU, let us consider the equation 

{^zUG + imz))=z\wk), (15) 

on Hg- We find that 

oo 
oo 

= ^z"t/aii,o)fc, (16) 

n=l 

which converges absolutely and uniformly on the disc \z\ < 1, and hence |$(2^)) 
is analytic in the disc \z\ < 1. 

Let Tio be the subspace of Hg spanned by the eigenvectors of U that are 
contained in Hg, i-e. eigenstates of U that have their support in the graph G. 
These are the analogs of bound states in conventional potential scattering, and 
in that case, the bound states are orthogonal to the scattering states. A similar 
situation obtains here. Let Hi be the orthogonal complement of Ho in Hg, let 
Pi be the orthogonal projection onto Hi, and let Ui — PiU . It is easily seen 
that \wk) and all of its images under f/i, are orthogonal to Hq. We now have 
that |$(2:)) can be expressed as 

oc 

|$(z)) =^z"+i[/,"|u;fc), (17) 

n=0 

takes its values in Hi, is unique, and is analytic for \z\ < 1. The solution can 
be extended beyond the disc \z\ < 1 by making use of the fact that the operator 
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[lui — zUi) on Til has no eigenstates in the closed disc \z\ < 1. Because this is 
an operator on a finite dimensional space, it has a finite number of eigenvalues, 
one of which is closest to z = 0. Let the magnitude of this eigenvalue be r. 
For \z\ < r, the inverse of {I-Hi — zUi) on Tii exists and is analytic [TS], and 
therefore, for \z\ < r we have that \^{z)) exists and is analytic. 

Summarizing, we have that Eq. P3|) has a unique solution, = |V^G,fc): 

with values in Tii, which is analytic on a domain \z\ < 1 + e, for some e > 0. 
Furthermore, we have that 

tf\z) = ,(0,l|C/(|V'G,fc(^)) + |l,0)fe), (18) 

so that t^l'\z) is analytic for \z\ < 1 + e, because \'4>G,k{z)) is. We also note that 



3 Spectral results 

Let 7ik be the closed, linear, [/-invariant subspace of 7Y = L'^{E) spanned by 
the states U^\j,j — for j = 1, 2 . . . and I an integer. This is just the space of 
states generated by incoming states on the fc*^ tail. We want to construct the 
spectral representation of U on Tik- We have that 

1. U\Me''))=e-''\Me'')) 

2. = 1 

3. {'ijjk{e''^)\v) = for any |w) G Tio ■ 

If \v) e L^{E), then {ipk{e^^)\v) is a continuous function of 9, so 

^ r d9\Me"')){Me''')\v) e L°-{E). (19) 
Zn 







The arguments in then show that for any \y) G Hk and for / any complex- 
valued continuous function on the unit circle (which we shall denote by C), 
that 

{y\fiu)\y) = ^ r defie-^'mkiOMl'- (20) 



2^ Jo 

Therefore, by the Riesz-Markov theorem, there exists a unique measure, /iy, a 
Borel measure on C, called the spectral measure associated to \y), such that 

{y\fiU)\y)^ r dt.y{9)f{e~^% (21) 
Jo 

and hence 

, ^ \iMO)\y)\\, 

27r 
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Now define the operator Vk : Hk ^^(C) by 



Vk\v) = -^me'')\v). (23) 



We then have the foUowing theorem [M]: Vk is a unitary operator from Hk to 
Li^iC) such that for any g{d) G L^{C), we have that 

VkUV,-'g{e) = e-'^g(0). (24) 

Hence, we have constructed part of the spectral decomposition of U. The U- 
invariant subspaces Tik are orthogonal for different values of fc, and they are 
all orthogonal to Tio- We have, in fact, that L'^{E) ^ Ho ®k=i ^fc EH- This, 
then, completes the spectral decomposition of U. 



4 Properties of the transmission coefficients 

The transmission coefficients, t^p {9) describe the behavior of a particle that 
starts on the k^^^ tail and scatters into the tail. Many properties of a quantum 
walk can be found directly from these functions. 

Suppose we start a walk in the state |0, l)k, and after each time step we 
measure the edge |0, l)j in order to see if the particle has arrived there. The 
probability that we find the particle there after n steps, but did not find it there 
for any of the previous n — 1 steps, which we shall denote by qf\n), is [14] 



(n) = U(0,l|[/"|l,0), 



(25) 



This probability can be expressed in terms of tj''\0) as follows. We have al- 
ready seen that \ipkiz)) is analytic in a region including the unit disc, and an 
examination of Eq. (jlip shows that it vanishes when z = 0. This implies that 
the functions t^^^ [z) and 'ipG.kiz)) are also analytic and vanish at z = 0. This 
implies that 



|l,0)fc 



1 



dz-\iJkiz)) 



1 

2^ 



de\MS)). 



(26) 



We then have 



so that 



fe(0,l|C/"|l,0), 



q) '(n) 



1 

2^ 



1 

2^ 



dOe' 



lOAk) 



(0), 



d0e- 



10, (k) 



(27) 



(28) 



This also gives us that the probability to find the particle on the j*"^ tail at some 
step is given by 



p(fe) 

" j,out 



oo 



(fe) 



in) 



1 

2^ 



de\tf\9)\ 



(29) 
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This probability, and other probabilities of interest, can also be expressed in 
terms of tj''\z) by means of contour integrals [T3]. 

The transmission coefficients also satisfy an orthogonality relation. In order 
to derive it, first define 

\MO)) = \i^G,km + \hO)k. (30) 

We then have that 

(Mmkie)) = {iJGA0MGMm + (31) 

If we apply U to 10^(0)), wc find 

K 

U\Mff)) = e-^'UoMO)) + E 4'^ (^)IO' 1)^-]' (32) 
from which it follows that 

K 

{Md)\MO)) = i^cimii'GMe)) + T.*?mf\o). (33) 

Comparing Eqs. (PTjl and ([55)1 . we see that 

K 

Y^tf{ertf\e) = 5u,i. (34) 

Let T{dG) = (|0, 1),,1 <l < K) and Lo{dG) = (|l,0)fe,l < fc < i-T) be the 
edges pointing out of G and the edges pointing into G respectively. Let T{dG) 
and ri(9G) be the corresponding subspaces of spanned by these edges. If we 
now define our scattering matrices 

S{e) : n{dG) T{dG) (35) 

by 

K 

S{\l,Q)u)^Y.^f\e)\Q,l)u (36) 
(=1 

our calculation shows that for each value of 9, S{9) is an isometry. We also see 
that 

\tf{9)\'de (37) 

is equal to the probability that a particle which starts at |1, 0)^ exits G at |0, l)j 
into the j*'' outgoing tail Yj. 

We will now investigate how this scattering matrix behaves when we subject 
r to an automorphism, and how it behaves under reversal. Let $ ~ (/, F) 
be quantum automorphism on T which induces a permutation tt^j on the edges 
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in u){dG) and a permutation tt^ on the edges of T{dG)). These permutations 
actually permute the corresponding tails. It is easy to see that 

So'Ka,=irroS, (38) 

if we extend tTi^ and tTt to be linear isometrics on n{dG) and T{dG) respectively. 
If we write tTojII, 0)fc = |l,0)7r^(fe) and tTt-IO, l)j = |0, l)^^(j) we see that 

The discussion of the effect of reversing the graph is somewhat more com- 
plicated. First we have to extend the reversing construction to graphs with a 
free quantum structure. Let F be a given Eulerian graph with incoming (out- 
going) tails Xi, ...,Xji (Yi, ...Yk) attached to G at vertices vi, ...vk [ui, ...uk) 
respectively. Let F^ be the reverse of F. We essentially reverse the orienta- 
tion on all the oriented edges on F which means that we swap the incoming 
edges for outgoing ones and vice versa. More specifically let R be the revers- 
ing map we defined in section 2 which reverses the orientation of each edge. 
Then i^f ) is the fc*^ (i*^) outgoing (incoming) tail of F^, which has the 
same vertices as Xk (Yj) but with the orientation on the edges reversed. If 
(/fe), where I = 1,2,... {{rrij), where m = 1,2...) are the vertices Xk (Yj) then 

the edges of Y^'^^ (X^) are |(Ofc,(Z + l)k) = + 1)^^ = R'Hll + hl)k) 
{\{m+l)j,{m)j} = |rn+ l,m)5^^ = R-Wm,m + 1)^). Let = R-^U'^R 
be the induced quantum structure on F^^^ Then 

c/(«)|z - 1,0^^ = R-'u-'\i,i- i)k = R-'\i + lj)k = \l,l + l)f ^ (40) 

Similarly we can show that 

U^^^\m+l,m)f^ = \m,m-l)fK (41) 

Therefore the induced quantum structure on F^^^ is free. 

Let Hg and 7i[^'' be the Hilbert spaces generated by the interior edges of 
G and G'-^^ respectively . Let Pq and Pq^' be the orthogonal projections onto 
Tic and T-Lq , respectively. It is easy to see that 

PgoR^Ro P^^ R-^oPg = P^^ oR-\ (42) 

Let Uq^^ = Pq^^ o {/(^) . We will now explicitly construct 

IV-f ^(^)) = W^-\z)) + |Vg](^)) + E(t(^))P^(^)k^f (^)), (43) 

1=1 

oo 

IV'SW) = E^"(4''V|l,0)f^ (44) 



where 
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formula (18) shows that 

= ((0,l)fV^''^[l^g]W> + |l,0)f^]) 

= (i?(0,l)f'|i?t/(«)[|V^g](z)) + |l,0)f)])* 

= (C/(l,0)fc|i?[|V^g](z))] + |0,l),r. (45) 



Now 

(C/(l,0),|i?(z"(4^))"(l,0)f)) 
This impUes 

which gives us the desired formula 



= (z"(C/Pg)"C/(1,0),|(0,1),) 

= (z"C/(C/G)"(0,l)fe|(0,l),). (46) 

= ([/^G,fc(^)|(0,l),), (47) 



(<(^))i^')(z) = (C/(|^G,,(^) + (1,0),)|(0, l),r = if (48) 

This enables us to compute the transmission amplitudes for T'^-^^ from those of 
our original graph F. If S{9) is the scattering matrix for F then S*''''{9) is the 
scattering matrix S'(^)(6') for F(«). Further more we see that 

K 

i?o^(«)(0)oi?-i|O,l)fc = ^4')(0)*|O,l)i (49) 

1=1 

It is because of Eq. ([T7|) that we need to choose R to be conjugate linear. 
The right hand side is clearly a conjugate analytic function of z, so the fact that 
R interchanges the complex structures forces the left hand side to be conjugate 
analytic too. When we defined the reverse graph Gr, we defined the mapping R : 
L'^{E}i) i^(i?) to be a conjugate linear isometry. If instead we replace L^(i?i?) 
withL2(^;^)^^ the space with the conjugate complex structure and Hermitian 
inner product, R becomes a complex linear isometry from L^{Eii)c onto LF'{E). 
If we follow the argument above with these changes we get the formula 

{t^^))^i^{z)=tf\z)\ (50) 

which seems to contradict analyticity until we recall that we are dealing with 
L^{Efi)c- So the "z" on the left hand side refers to the complex structure which 
is the conjugate of the original one on L?'{E). Hence, the more functorial formula 
above is as it should be. If we followed through in this vein we would finally 
arrive at Eq. (|48|) . We chose to avoid this additional complication because our 
concern was to give formulas to compute the transition amplitudes. 

5 Add a handle 

We now want to begin our study of how changing a graph modifies its transmis- 
sion amplitudes. Let us first consider what happens when we replace a pair of 
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tails, one incoming and one outgoing with a single edge going from the attaching 
vertex of the outgoing edge to the attaching vertex of the incoming edge. This 
new graph V will still be Eulerian. It will also inherit a free quantum structure 
U' from the original free quantum structure [/ on F. 

Let us be more specific. Let F = (G, {vi, Vk), (ui, uk)) be an Eulerian 
graph with tails. Let Xi, ...Xk be the incoming tails attached at the vertices 
vit..,vk respectively and let Yi^...,Yk be the outgoing tails attached at the 
vertices u\, uk respectively. Let F' = (G", (112, ■■■vk), (^2, uk)) be the new 
Eulerian graph with tails. The graph G' has the same vertices as G. Its edges 
are those of G plus a new one |w,i,i;i), which is an oriented edge from ui to vi. 
Its incoming (outgoing) tails are X2, ...Xk iY^, ...^Yk) with the same attaching 
vertices as in F. We see that = r„ and u>'^ = uiy for vertices v not equal to 
Ml or Vi. The prime refers to F'. The edges of uj[.^ are the same as those of 
LOvi with 1 1,0)1 = |li,i'i) replaced by \ui,vi). The edges r^^ are the same as 
those of with |0, l)i = |ui,li) replaced by \ui,vi). We can now see how 
the original free quantum structure C/ on F induces a new one, U', on F'. If a 
vertex v is unequal to ui or vi then = Uy. If the vertex is vi or ui this is 
still the rule if we consistently replace |l,0)i in uiy^ with \ui,vi) and |0, l)i in 

with \ui,Vi). So 

(7:Jmi,i;i) = (7„|L0)i, (51) 

which is unambiguous unless ui = Vi. In that case Uy\l,0)i may contain a term 
proportional to \ui, (l)i). In defining [/^^ \ui,vi), we simply take the expression 
for ?7-„|l,0)i and replace \ui, (l)i) by |?ii,?;i) wherever it occurs. 

Our aim, now, is to show how to simply compute the transmission amplitudes 
rj*^^ (z) of the new configuration F' from the transmission amplitudes tj'^ {z) of 
the original F. We will adopt the notation of section 2. Let 

mz)) = Wk-{z)) + li^G'Az)) + Sf=2rf (^)k,.+ (^)), (52) 

be the generalized eigenstate of U' which satisfies the equation 

zU'l^iiz)) = \i,',{z)). (53) 

We want to represent |V'G',fe(^)) in the form 

li^G'Az)) = IV'G,fc(^)) + (ik{z)i\ui, vi) + |Vg,i(^))), (54) 

where \ipG,k{z)) is that part of the generahzed eigenstate of F that is supported 
in G, and au (z) is a function to be determined. We have that 

K 

zU'{\u„v,) + |^G,i(^))) = |^G,i(^)) + 4'^(^)|Mi,^;i) + ^t«(z)|0, 1)„ (55) 
and 

K 

zU'{\l, 0)k + |VG',fc(^)) = |V'G',fc(^)) + Yl ^i'^ (^)IO' (56) 

j=2 
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where, as previously mentioned, the rj'^-' (z) denote the transmission amplitudes 
for r'. If we represent |i/'G',fe(-z)) as in Eq. (f54l we see that 

zU'{\l,0)k + HgA^)) + ak{z)i\u,,v,) + HjgAz)))) 

K 
2 

K 

+t'^l\z)\u,,v,) + Y,tf\z)%l),). (57) 
i=2 



Comparing this result with Eq. (|56p . we can see that 

r^{z) = tf\z) + ak{z)tf\z) 

akiz) = t^^\z)+akiz)t[^\z), (58) 



so that, finally, 

rfH^)=tf(z) + i^^fl||^^, (59) 

^ ^ i-4^(z) 

all of which hold in a neighborhood of the unit disk by the discussion in section 
2. 

We can iterate this procedure by connecting an outgoing tail to an incoming 
tail as we did above, one pair at a time, thereby adding several "handles". We 
can also use this method to accomplish this in a single step in the following way. 
Let us try to add L handles by splicing Yj to Xj 1 < j < L respectively and 
form new edges (handles), ei = (ui^vi) 62 — (^2,^2) ...cl = {ul,vl). The new 
generalized eigenstate coming from an incoming wave traveling along one of the 
remaining incoming tails is 

K 

W,{z))^a^-{z)^\4>G',k{^))+ E rf\z)\a,.{z)), (60) 

j=L + l 

where 4'G'.k{z)) is of the form 

L 

li'G'Az)) = |i/'G,fc(^))+IZaM(^)[k;,«i) + lfe(^))], (61) 
1=1 

where the functions ak.i{z) are to be determined. Exactly as before, we have 
that 

K 

zC/'(|l,0)fe + |^G',fc(^))) = |V'G',fe(^))+ (62) 

j = L + l 

where rj'^^ (z) denotes the appropriate transmission amplitude for T' the new 
Eulerian graph with tails where we replaced the first L pairs of tails with the 
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corresponding handles. As before 



1=1 



zU'{\l,0)k + |^G,fc(^)) + ^afe,i[|l,0)i + IVg/^))]) 

1=1 

j=L+l j=l 
K L 



(63) 



j=L+l 3 = 1 

These equations lead to two sets of equations. The first set, L in number, is 

L 



(64) 



where 1 <l < L. This set of equations can be solved by Cramer's rule to give 

expressions for a^j as rational functions of (z) and (z) in some neighbor- 
hood of the origin in the complex plane, because, the transmission amplitudes 
of r vanish at the origin. The second set of equations is 



(65) 



1=1 



iov L + 1 < j < K. If we substitute the solutions to the first set of equations 
into the second set we get our desired result, which expresses the transition 
amplitudes of V as rational functions of the transition amplitudes of V. The 
formula holds on a neighborhood of the unit disk in the complex plane. The 
rational functions themselves only depend upon the numbers K and L. 



6 Cut a handle 

Let us start with F an Eulerian graph with tails. Let X2-,...-,Xk {Y2,...,Yk) K—1 

incoming (outgoing) attached at the vertices V2,...,vk {u2t..,uk) respectively. 
Let e = (ui, vi) be an oriented edge in G. We wish to replace the edge e with a 
pair of tails, an incoming tail Xi attached at vi and an outgoing tail Yl attached 
at ui. The new vertices on Xi the incoming edge will be denoted 1^:, 2a;,... and 
the corresponding edges will be denoted |fc, k — l)x = \kx, {k — l)a;). The new 
vertices on Y\ the outgoing edge will be 1^,2^... and the corresponding edges will 
be k + l)y — |fcy, [k + l)j,) fc = 1, 2, .... We will sometimes denote ui as Oj, 
and denote V\ as 0^ respectively. Let G' be the graph with the same vertices 
as G and the same edges as G with the edge {u\,v\) removed. The new graph 
with tails T' ~ (G' ,{vi,V2, ■■■,vk),{ui,U2, ...tUk)) has all the vertices of T as 
well as the new vertices kx and ky for the two new tails. The edges of F' are 
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those of r with {ui,vi) removed and the tail edges \k,k + 1)^ and \k + 
k = 1,2, ... added. The original free quantum strueture [/ on F induees a new 
free quantum structure U' on V as follows. First Uj, \k — 1, k)y = \k,k + l)y 
and [/^ \k + l,k)x = \k,k — 1)^. must hold in order that U' bo free. If v is a 
vertex of G' which is not equal to ui or vi then U!^, = f/„ . In order to deal with 
vi and Ml we note as in the last section that the edges of lo'^-^ (f^J are the same 
as those of lu^-^ {tui) with the edge (ui,fi) of T replaced by |l,0)a; = |la;,wi) 
(|0, l)y = \ui, ly). Thus the rule = Uy still holds for vertices vi and ui if we 
consistently replace \ui,vi) in coy^ {tui) by |l,0)a; (|0, 1)^) so 

U:,^\l,0), = U,,\ui,vi). (66) 

This is unambiguous unless Ui = Vi where we replace \ui,Vi) by |0, 1)^^ in the 
formula. 

This operation is clearly the inverse of our add a handle procedure. If we 
create a pair of tails Xi and Yi from an edge {ui,vi) as above and the splice 
the new tails together to reform the edge {ui,vi) according to our add a handle 
prescription wc; end up with the same free quantum structure on the same graph 
with tails. Similarly if we start by first adding a handle by splicing a pair of 
tails, and then cut this new handle by the rules above we again arrive back 
where we started from. 

Let t^j'^\z), 2 < j, k < K,he the transmission amplitudes for T, our original 

graph with tails, and let T^''\z), 1 < k,j < K, be the transmission amplitudes 
for the new configuration T'. We wish to compute the functions t'^'^\z) from 

the attributes of F, which include the functions t'j^\z), as simply as possible. 

We begin by applying the results of the previous section to the graph F', as F 
is obtained from F' by adding a handle. This immediately gives us that 



tf\z)=T^'\z) + ^' (67) 



for 2<j,k<K. If we could find the functions Tj'^\z) and tI^\z), for 1 < i < 

K, then we could solve these equations for the remaining Tj''\z). Note that 
these are the transmission amplitudes associated with the new tails, Xi and Yi . 

In order to describe how to calculate these transmission amplitudes from F, 
we define the generalized eigenstate on F' 

mz)) = \a,-{z)) + \i>G',i)+Y.Ti'\z)\aj,{z)), (68) 

i=i 

This is the formula for the generalized eigenstate generated by an incoming wave 
along Xi, where \ipG',i) is the part of this eigenstate that is supported on He- 

Let us note that Tic is the direct sum of Tic and the one-dimensional subspace 
consisting of multiples of \ui,vi). Let Pc denote the orthogonal projection 
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operator onto Wg'j and let Uc = Pg'U. Considered as a state on F', IV'ca) 
satisfies the equation 

zPg'U'H^G'Az)) + |1,0),) - \^g'a{z)), (69) 
and this implies that considered as a state on T it satisfies 

zPg'UHiPg'Az)) + \ui,vi).,) = IVgmW). (70) 
The solution to this equation is 

oo 

|V'gm(^)) = E^"4'Vi,^^i)- (71) 

Tl=l 

Note that this equation contains only quantities defined on the original graph, 
r, so that it implies that |'i/'G',i(^)) can be calculated from the initial graph. 
Once we have found \ipG',i{z)) we can substitute it into the equation 

7;W(^)^, (0,1||[/'(^gm(z) + |1,0).), (72) 

defined on T' , to find T^^\z). 

Our remaining task is to find T^''\z), for 2 < j < K. This can be done 
by looking with the reverse graph and its induced quantum structure and 
following a procedure analogous to the one we followed for T. We first find 

oc 

\^B',iiz)) = ^z"(P«C/(^))«|x;i,«i), (73) 

n=l 

where Pq, projects onto the subspace spanned by the states corresponding to 
all of the edges of G^^^ except \vi,ui). Once we have \iIjq, i{z)), we can use it 

to find Tj^'''"^\z) in the same we we used \ipG',iiz)) to find Tj^\z). Finally, we 
can make use of the relation 

7;(«)W(,)=r(-'")(z), (74) 

to find T^-'\z). This means that all of the quantities in Eq. (|^ except t'-^\z) 
for 2 < k < K are then, in principle, known so that they can be used to find 
the transmission amplitudes for the cut graph. 

7 The splice 

We can use the results of section 5 to find the transmission amplitudes of a graph 
that is the result of splicing two other graphs together. Let us assume that the 
graph G, with entering tails Xi, . . . Xk attached to G at the vertices vi, . . .vk, 
and exiting tails Yi, . . . Yk attached to G at the vertices mi, . . . uk, can be broken 
into two disjoint pieces, Gi and G2. Gi contains the vertices to 
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which the tails Xi, . . . X^^i are attached and the vertices ui . . . tiL-i to which 
the tails Yi. . . Y^-i are attached. Similarly, G2 contains the vertices Vl, ■ ■ ■ Vk, 
to which the tails Xl, ■ ■ ■ Xk are attached and the vertices ul ■ ■ ■ uk to which 
the tails Yl,...Yk arc attached. We now remove the tails Yi and Xl and 
replace them by one directed edge from ui to vl- This is simply a specific 
example of the adding a handle construction. 

Let us now apply the results of the previous section. We first note that in 
the original graph, a particle entering Gi would never exit from G2, so tj'^^ = 
if Vk S Gi and uj G G^ ■ Similarly, a particle entering G2 would never exit from 
Gi, which implies that t^'^^ = if -Ufc € G2 and uj G Gi. Making use of this 
result, we have from section 5 that 



T, 



^ 1^ tj'\z) otherwise 

Therefore, we can use the add-a-handle construction to construct a quantum 
walk on a larger graph from walks on smaller ones. 



8 Comparing graphs 

We can use the techniques developed here to compare two graphs by doing 
a kind of interferometry with them. Suppose that Gi has only two tails, an 
incoming tail Xi attached to vertex vi and an outgoing tail Yi attached to ui . 
Similarly, G2 also has only two tails, an incoming tail X2 attached to V2 and an 
outgoing tail Y2 attached to U2- We would like to determine whether Gi and 
G2 are the same or different, with the additional constraint that if they are the 
same, vi in Gi should correspond to V2 in G2, and Ui in Gi should correspond 
to U2 in G2. 

One way of attacking this problem is to put the two graphs into an arrange- 
ment like that shown in Figure 1. First, we remove the tails from the graphs. 
We now consider a vertex. A, with two incoming and two outgoing edges. The 
incoming edges are the initial edges of two incoming tails, Xia and X2a- One 
outgoing edge is attached to vi and the other is attached to V2. Similarly, we 
consider a second vertex, B, with two incoming and two outgoing edges. The 
outgoing edges are the initial edges of two outgoing tails, Yis and 125, and 
one incoming edge is attached to ui and the other is attached to U2. We shall 
denote by F be the graph with tails that we get when we connect Gi and G2 to 
the vertices A and B with the tails Xia, X2A, Yrs and Y2B, and by U the free 
quantum structure on F. 

What we have done is to create a larger graph from Gi and G2 in which the 
two graphs are in parallel. We shall choose the unitary operators corresponding 
to the vertices A and B in such a way that if we start the walk on the incoming 
tail XiA, and if the graphs Gi and G2 are identical, with vi in Gi corresponding 
to V2 in G2, and Ui in Gi corresponding to U2 in G2, then the transmission 
amplitude corresponding to the outgoing tail F2B will be zero. Therefore, if we 
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Figure 1: An arrangement for comparing the graphs Gi and G2. A particle 
starting a walk on the tail Xi will never exit on the tail Y2 if the graphs Gi and 
G2 are identical. 

start a walk on the tail Xia, and eventually find the particle on the tail Y2B, 
we can conclude that the two graphs are not identical. 

We now want to find the transmission amplitudes of the combined graph in 
terms of the transmission amplitudes of Gi and G2. In order to do this, we need 
to specify what happens at the vertices A and B. The states corresponding to 
the edges entering vertex A arc which is the initial edge of Xia, and 

|l2,^), which is the initial edge of X2A- The states corresponding to the edges 
leaving vertex A are \A, vi) and \A, ^2). The local isometry at vertex A mapping 
incoming to outgoing states has the following action 

111, A) ^ l=(\A,v,) + \A,V2)) 

|l2,A) ^ -^i\A,v^)-\A,V2)). (76) 

The local isometry at vertex B is essentially identical. The states corresponding 
to the edges entering vertex B are \ui, B) and \u2, B). The states corresponding 
to the edges leaving vertex B are |_B, li), which is the initial edge of Yis, and 
\B,V2), which is the initial edge of ¥23- The local isometry at vertex B mapping 
incoming to outgoing states is given by 

\ui,B) 

\U2,B} 



V2 
1 

7i 



(|S,li) + |S,l2)) 

{\B,U)-\B,l2)). 



(77) 
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The generalized eigenstate corresponding to the entering particle being on 
the tail Xia is 

2 

\Mz)) = ki-W) + IV'a,b(^)) +J2^^'\z)\aj+{z)), (78) 

where \iPa,b{z)) is the part between vertices A and B. The transmission ampli- 
tude t[^^ {z) corresponds to the tail Yib, and the transmission amplitude T2^^ (z) 
corresponds to the tail I2B • We would like to find these transmission amplitudes 
in terms of the transmission amplitudes for Gi and G2, which we shall denote 
by t'"il{z) and 4^2 (•^)' respectively. We can find the transmission amplitudes for 
r, and the state \iPa,b{z)) by making use of the equations 

2 

zU{\U,A) + \i,AM^))) = \i>AM^))+Y^r^'\z)\B,lj) 

zU{\A,Vk) + \Mz))) = \Mz))+t''^,l{z)\uk,B), (79) 

where fc = 1,2 in the second equation. The states Itpkiz)), for k = 1,2 are 
the internal parts of the generalized eigenfunctions of the graphs Gi and G2, 
respectively. We find that 

2 

rW(.) = il(tW(.)-iW(.)). (80) 

We can now clearly see that if the graphs Gi and G2 are the same, which implies 
that t^iliz) = t^ili^)' then T2'\z) = 0. Therefore, in this case the particle will 
never enter the tail ¥23- 

We have shown that we can compute the transmission amplitudes for the 
combined graph, F, in terms of those of Gi and G2, and that one of these 
amplitudes vanishes if Gi and G2 arc the same. Therefore, a quantum walk on 
the combined graph can be used to compare Gi and G2. An important issue is 
how efficient this procedure is. In particular, we would like to know how many 
steps the walk needs to make to determine whether the graphs are the same or 
different, and how the number of steps is related to the number of edges in Gi 
and G2. This will remain for future work. 

9 Conclusion 

We have presented a formalism that describes quantum walks on directed, Eu- 
lerian graphs. These graphs have incoming and outgoing tails, and we define 
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a free quantum structure on them that determines the motion of a particle on 
the graph. We found that the transmission amphtudes of such a graph are very 
useful in determining the behavior of a quantum walk on that graph. We have 
shown how the transmission amplitudes of altered graphs, i.e. graphs that are 
formed from an original one by adding or cutting an edge, can be found from 
the transmission amplitudes of the original graph. This allowed us, in some 
cases to find the transmission amplitudes of a graph in terms of those of its 
subgraphs. Finally, we showed how these constructions can be used to compare 
two graphs by constructing a larger graph from two smaller ones, which acts 
as a kind of interferometer. If the two smaller graphs are identical, a particle 
making a quantum walk on the larger graph can only emerge onto one of the 
two exit tails, but not onto the other. 
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